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ABSTRACT 

We have constructed a sample of radio-loud objects with optical spectroscopy from the Galaxy 
and Mass Assembly (GAMA) project over the Herschel-KYLAS Phase 1 fields. Classifying 
the radio sources in terms of their optical spectra, we find that strong-emission-line sources 
('high-excitation radio galaxies') have, on average, a factor ~ 4 higher 250-/im Herschel lu- 
minosity than weak-line (iow-excitation') radio galaxies and are also more luminous than 
magnitude-matched radio-quiet galaxies at the same redshift. Using all five H- ATLAS bands, 
we show that this difference in luminosity between the emission-line classes arises mostly 
from a difference in the average dust temperature; strong-emission-line sources tend to have 
comparable dust masses to, but higher dust temperatures than, radio galaxies with weak emis- 
sion lines. We interpret this as showing that radio galaxies with strong nuclear emission lines 
are much more likely to be associated with star formation in their host galaxy, although there 
is certainly not a one-to-one relationship between star formation and strong-line AGN activ- 
ity. The strong-line sources are estimated to have star-formation rates at least a factor 3-4 
higher than those in the weak-line objects. Our conclusion is consistent with earlier work, 
generally carried out using much smaller samples, and reinforces the general picture of high- 
excitation radio galaxies as being located in lower-mass, less evolved host galaxies than their 
low-excitation counterparts. 

Key words: galaxies: active - radio continuum: galaxies - infrared: galaxies 
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1 INTRODUCTION 

The relationship between AGN activity and star formation is a com- 
plex one. In order to maintain the observed black hole mass/bulge 
mass relationship, bla ck holes must grow as new stars form (e.g. 
iMagorrian et alj Tl998b and black hole growth should result in AGN 
activity. The generally accepted picture is one in which mergers 
trigger both AGN activity an d star formation (e.g. iGranato et all 
|2004 |Pi Matteo et all 120051) and in which the AGN activity, at 
some point, shuts down star formation by on e of a range of pro - 
cesses generally referred to as 'feedback' (e.g. lCroton et alj,l2006h . 
The microphysics of this process presumably invo lves the drivin 
of ou tflows either by lumin ous quasar activity (e.g . lMaiolino et al 
|2012|) and/or radio jets (e.g. lHardcastle et alll2012l) ; an understand- 
ing of how feedback operates in populations of galaxies is crucial 
to models of galaxy and black hole evolution. 

Radio-loud active galaxies form a particularly interesting sub- 
population of AGN in the context of this question. Firstly, they tend 
to reside in massive elliptical galaxies, traditionally thought to be 
'red and dead' with little or no recent star formation; secondly, 
the large amount of kinetic energy that they inject into their en- 
vironment means that they must both influence and be influenced 
by the galactic environment in which they are em bedded. There is 
in fac t long-standing observational evidence (e.g. iHeckman et all 
Il986t) that some powerful radio galaxies have peculiar optical mor- 
phologies, plausibly the results of mergers with gas-rich galaxies. 
In these systems, we might expect AGN activity and star forma- 
tion to go hand in hand, although the different timescales for star 
formation and AGN triggering will mean that they will not always 
be observed together; in radio-quiet systems, there may be several 
hundred Myr of delay betw een the starburst and the peak of AGN 
accretion dWild et allkoiOl) . In contrast to these objects, we know 
that other radio galaxies, often equally powerful when their kinetic 
powers can be computed, reside in the centres of rich cluster or 
group environments where, on the one hand, gas-rich mergers must 
be very rare, and, on the other, the duty cycle of AGN activity must 
approach 100 per cent to account for the nearly universa l detec - 
tion of radio sources in these systems (e.g. lEilek & Owenl l2006t) . 
In these objects, we would be surprised to see evidence for a direct 

link between AGN activity and star formation. 

It may be possible, as originally suggested bv lHeckman et alj 

dl986h , to understand these apparently contradictory results 
in the context of a two-population model of the AGN activity 
in radio galaxies. The two populations in qu estion probably 
corres pond quite closely to classes A and B of llline & Longain 
dl979h . now known as high-excitation and low-e xcitation ra- 
dio g alaxies (hereafter HERGs a nd LERGs: e.g. lLaing et all 



1 1994 ; Ijackson & Rawlirigsl 1 19971) . In recent years it has be- 
come clear that the differences between these objects are not 
simply a matter of emi ssion-line strength but extend to optical 
1 Chiaberge et"ail l2002n . X-ray dHardcast le. Evans, & Croston, 
20061) and mid-IR dOgle. Whvsong. & Antonuccil 12006 



Hardcastle. Evans. & Crostonl l2009Th In the vast majority of 



LERGiS there is no evidence for any radiatively efficient AGN 
activity, setting aside non-thermal emission associated with the 

1 The optical emission-line class doe s not correspond compl etely reliably 
to other indicators of AGN activity; see Hardcastle et al. ( 2009 ) for a discus- 
sion of some anomalous or intermediate objects and lRamos Almeida et alj 
b'l 1,1 ) for a particularly well-documented 'LERG' with a clear heavily 
absorbed, luminous hidden AGN. Emission-line classification clearly does 
not have a one-to-one relationship to radiative efficiency, but, for simplicity, 



nuclear jet (e.g. lHardcastle et alj 120091 and references therein); 
the AGN power output is primarily kinetic and we observe it 
only through the radiation of the jet and lobes and through the 
work they do on the medium in which they are embedded. On the 
other hand, the HERGs, which include the traditional classes of 
narrow-line radio galaxies with spectra like those of Seyfert 2s, 
the broad-line radio galaxies and the radio-loud quasars, behave 
like textbook AGN with the addition of jets and lobes. Although 
LERGs are more prevalent at low radio powers and HERGs at 
high powers, both classes are found across the vast majority of the 
radio power range and, where they overlap, there is often no way 
of distinguishing between the radio structures that they produce. 

The reason for the fundamental differences between the AGN 
activity in these two classes of radio source is not clear, but one 
proposal is that the difference s arise because of diffe rent fuelling 
mechanisms. In this scenario dHardcastle e7aill2007l) the LERGs 
are fuelled directly from the hot gas halos of their host ellipticals 
and the groups and clusters in which they lie, while the HERGs 
are fuelled, often at a higher rate, by cold gas, presumably brought 
into the host elliptical by mergers or interactions with gas-rich sys- 
temfl Because the LERGs dominate the population at low power 
and low redshift, this allows a picture in which nearby radio-loud 
AGN are driven by accretion of the hot phase and are res ponsi- 
ble for balancing its radiative cooling (e.g. lBest et al| [2006) while 
still allowing for merger- and interaction-driven radio-loud AGN 
at higher radio luminosity and/or redshift. This model makes a 
number of testable predictions. LERGs will tend to be associ- 
ated with the most massive systems, will therefore tend to inhabit 
rich environments, and will largely have old stellar populations; as 
a population, they will evolve relatively slowly. HERGs can oc- 
cur in lower-mass galaxies with lower-mass black holes, provided 
that there is a supply of (cold) fuel: we therefore expect them 
to be in less dense environments, to be associated with merger 
and star-formation signatures, to be in less evolved, lower-mass 
galaxies and to evolve relatively fast with cosmic time (since the 
merger rate was higher in the past). Many of these predictions 
have been tested. There is some evidence, particularly at low red- 
shifts, for a difference in the enviro nments and the masses of the 
host g alaxies of LERGs and HERGs dHardcastlej|2004lTasse et all 
I2008L Ching et al., in prep.), and there is strong evidence, also 
at low redshifts, for differences in the host galax y colours in the 
sense expected from th e model described ab ove dSmolcidl2009l ; 
iBest & Heckmanl 1201 2| ; I Jans sen et all l2012h . There is strong ev- 
idence for an increased fraction of signatures of merger or inter- 
action in the galaxy morphologies of the H ERGs with respect to 
the LERGs dRamos Almeida et al , 201 lbl) and with respe ct to a 
background galaxy population 1 Ramos Almeida et alll2012h . And, 
most importantly from the point of view of the present paper, there 
is direct evidence for different star-formation histories in the hosts 
of HERGs and LERGs, in the sense predicted by the model, i.e. 
that HERGs show evidence for more recent star formation both at 



in this paper we will continue to refer to LERGs and HERGs as though they 
represent the archetypes of their population. 

2 It is not yet clear whether the difference in the AGN results from the dif- 
ference in the temperature of the accreted material, as proposed by Hardcas- 
tle et al., or simply from the lower accretion rates as a fraction of Eddington 
expected for massive black holes being fed at something approx i matin g the 
Bondi rate in the LERGs, as in the m odels of Merloni & Heinz 1 200§|) and 
as argued bv lBest&Heckmad 120121) ; Mingo et al. (in prep) will discuss 
this question in detail. However, the answer to this question makes very 
little difference to the predictions of the model. 
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low r edshift dBaldi & CapettHl2Q08l) and at z ■ 

lioToh . 

Studies of the star formation in the different classes of ra- 
dio galaxy have until recently been limited in size because of 
the techniques and samples used (e.g. HST imaging by Baldi & 
Capetti, analysis of optical spectroscopy by Herb ert et al.). Only re- 
cently have large samples b egun to be analysed dBest & Heckmanl 
[201 21; Uanssen et al. [ l2012h and so far this work has been based 
only on optical colours at low redshift. Mid-infrared observations 
with Spitzer provide some evidence that i ndividual HERGs may 
have strong star formation (e.g. Cygnus A. IPrivon et alj|2012l) but 
systematic studies of large samples have generally shown that 
the luminosity in the mid-IR is dominate d by emission from the 
AGN itself, by way of the dusty torus (e.g. lHardcastle et"aill2009l ; 
iDicken et al. , 20091); det ailed mid-IR spectroscopy in small sam- 
ples 1 Dicken et all 1201 2|) has shown that there is not a one-to-one 
association between star formation signatures and AGN activity, 
but this type of work cannot easily be extended to very large sam- 
ples. However, observations of cool dust in the far infrared (FIR) 
should, in principle, provide a very clear way of studying star for- 
mation, which should be uncontaminated by AGN activity, since 
the emission from the dusty torus of th e AGN is found to peak 
in the rest-frame mid-IR (e.g. lHaas et al.Ll2004l) . FIR observations 
can be carried out simply for large samples, and the method can 
extend to relatively high redshifts, with the only contaminant being 
emission from diffuse dust heated by the local interstellar radia- 
tion field rather than by young stars (at least until redshifts become 
so high that rest-frame mid-IR torus emission starts to appear in 
the observer-frame FIR bands). Earlier work on far-infrared/sub- 
mm studies of star formation in samples of radio galaxies neces- 
sarily concentrated on high-redshift objects, in which emission at 
long observed wavelengths (e.g. 850 flm, 1.2 mm) corresponds to 
rest-frame w avelengths around the expected peak of ther mal dust 
emission (e.g. lArchibald et alll200ll ; lReuland et a"Hl2004l) and thus 
applied only to very radio-luminous AGN. Much larger and more 
local samples can be stu died using the Herschel Space Observa- 
tory dPilbrattet all 1201 Cj|) and in particular by wide-field surveys 
such as the Herschel Astroph ysical Terahertz Large Area Survey 
(H- ATLAS ; lEales et all[2Qloh. 

In an earlier paper ( Hardcastl e et alll2010l . hereafter H10), we 
studied the FIR properties of radio-loud objects in the 14-square- 
degree field of the Science Demonstration Phase (SDP) dataset of 
H-ATLAS, and showed that, as a sample, their FIR properties were 
very similar to those of normal radio-quiet galaxies of similar mag- 
nitude; however, our sample size was small and we were not able to 
classify our radio-loud objects spectroscopically. The full 'Phase 1 ' 
ATLAS dataset, consisting of three large equatorial regions, gives a 
field almost twelve times larger (161 square degrees). Our work on 
radio galaxies in the Phase 1 dataset is divided between two papers. 
Virdee et al. (2012; hereafter V12) use the same sample selection 
process as H10, but use the much larger sample available from the 
Phase 1 datasets to investigate the relationship between radio galax- 
ies and normal galaxies in more detail, dividing the radio-loud sam- 
ple by properties such as host galaxy mass and radio source size. 
In the present paper, we select our sample so as to be able to clas- 
sify our radio sources spectroscopically, using data derived from 
the Galaxy and Mass Assembly project l lDriver et alj|2009l 1201 ll ; 
hereafter GAMA) and search for differences in the far-infrared and 
star-formation properties of HERGs and LERGs. 

Throughout the paper we use a concordance cosmology with 
H = 70 km s -1 Mpc~', fl m = 0.3 and £1 A = 0.7. Spectral index 
a is defined in the sense that S <x v~ a . 



0.5 faerbert et all 2 SAMPLE SELECTION AND MEASUREMENTS 



2.1 The GAMA sample 

The GAMA survey is a study of galaxy evolution using multiwave- 
length data and optical spectrosopic data. In phase I of GAMA, tar- 
get galaxies are drawn from the SDSSDR6 photometric catalogue 
in three individual 12° x 4° rectangles along the equatorial regions 
centred at around 9, 12 and 15 hours of right ascension. A r-band 
magnitude limit of 19.4 was used for the 9 and 15-h fields while 
the 12-h field had a deeper 19.8 mag limit jDriver et alll201ll) . The 
H-ATLAS Phase I data is taken from regions corresponding closely 
to these three fields. Reliable spectroscopic redshifts from previous 
surveys (e.g. SDSS, 6dF Galaxy Survey, etc.) were used for GAMA 
sources that had them. Those without reliable spectroscopic red- 
shifts from previous surveys were spectroscopically observed on 
the Anglo- Australian Telescope (AAT). 

We built a sample of candidate radio galaxies by cross- 
matching the Faint Images of the Radio Sky at Twenty-cm (FIRST, 
lBeckeretailll995h catalogue (16 July 2008) with optical sources 
(i < 20.5 mag, extinction corrected) fr om the Sloan Digital Sky 
Surve y Data Release 6 (SDSSDR6; lAdelman-McCarthv et all 
l2008h in all GAMA regions. The full details of the cross-matching 
will be described by Ching et al. (in prep.), but a short summary is 
provided here. The cross-matching firstly involved grouping FIRST 
components that were likely subcomponents of a single optical 
source (e.g. the core and lobes of a radio galaxy). The optical 
counterparts for the groups were matched automatically if they sat- 
isfied certain criteria based on symmetries of the radio sources, 
and/or manually when groups were more complex, by overlaying 
SPSS images with FIRST and NRAO VLA Sky Survey (NVSS; 
ICondon et all 1 19981) contours. Groups that appeared to be sepa- 
rate individual radio sources were split into appropriate subgroups 
matched to their individual optical counterpart. All FIRST com- 
ponents that were not identified as a possible subcomponent were 
cross-matched to the nearest SDSS optical counterpart with a max- 
imum separation of 2.5 arcsec. This process gave us a sample of 
3168 objects with radio/optical identifications. Some of these ob- 
jects, predominantly at low redshifts, had spectra from the SDSS 
spectroscopic observations; GAMA does not re-observe such ob- 
jects. Others were part of the GAMA main sample. To increase the 
spectroscopic sample size, we identified galaxies that were not part 
of the GAMA main sample (internal data management unit Tiling- 
Catvl6, SURVEY.CLASS^ 1), and observed some of them as 
spare-fibre targets during the main GAMA observing programme 
(see Ching et al., in prep., for more details). The resulting sam- 
ple, by construction, contained only sources with usable spectra 
and spectroscopically determined r edshifts (nQ > 3, f rom GAMA 
data management unit SpecCatv08; lDriver et alj201lh . and is flux- 
limited in the radio, with a lowest 1.4-GHz flux density around 0.5 
mjy and most sources having flux density above 1.0 mjy, as a result 
of the use of FIRST in constructing the sample. There were 2559 
sources with spectroscopic redshifts in this parent sample. 



2.2 Spectral classification 

Spectral classification of the objects with radio/optical identifica- 
tions was carried out by inspection of their spectra. A detailed de- 
scription of the process will be given by Ching et al. (in prep.); 
here we simply summarize the steps we followed. The emission 
line measurements used in this paper for GAMA spectra were made 
from the Gas and Absorption Line Fitting (GANDALF; ISarzi et alj 
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120061) code as part of the GAMA survey (see lHopkins et alj|2012l 
for a description of the GAMA spectroscopy and spectroscopic 
pipeline), while for SDSS spectra we used the measurements from 
the value-added MPA-JHU emission-line measurements derived 
from SDSS DR"0- 

Both of these measurements fit the underlying 
stellar population before making emission line measurements, and 
hence take into account any stellar absorption. Only high-quality 
GAMA spectra were used. 

We firstly removed Galactic sources by imposing a lower 
redshift limit of z > 0.002. Such objects are classified 'Star' and 
play no further part in the analysis in this paper. Next, we visu- 
ally selected objects with broad emission lines; these are classified 
'AeB' in this paper, and are broad-line radio galaxies or radio-loud 
quasars. 

Galaxies that are within z < 0.3 and have 1.4-GHz luminos- 
ity (hereafter Z-1,4) below 10 24 W Hz -1 have a high probabil- 
ity of having star-form ation dominated radio emission (see e.g. 
iMauch & Sadlerll2007l) . For all such objects having [Olll], [Nil], 
Ha and Hj6 emission lines detected with a sig nal-to-noise ra- 
tio > 3, we used a simple line diagnostic (BPT; iBaldwin et all 
198 lh to classify 'pur e star-forming galaxies' as classified by 
Kauffmann et al \ j2003h . These are classed as 'SF' in the following 
analysis. However, as pointed out bv lBest et all d2005h . line diag- 
nostics alone are not enough to ensure a clean sample of radio-loud 
AGN, since the emission lines may arise from a radio-quiet AGN, 
while the detected radio emission might arise from star formation in 
another region. In addition, the lines required for BPT analysis are 
not available at z > 0.3. We therefore also classified as 'SF' any ob- 
ject whose Ha and 1.4-GHz radio emission placed it within 3ct of 
the rel ation between these two quantities derived bv lHopkins et alj 
d2003h for star-forming objects. We emphasise that 'SF' objects are 
not discarded from the analysis at this stage - therefore nothing in 
this classification prejudices the results of the H-ATLAS analysis. 

Finally, we expected the remaining galaxies to be a reason- 
ably robust sample of radio-loud AGN, possibly contaminated by 
z > 0.3 and/or extremely luminous (L1.4 > 10 24 W Hz -1 ) star- 
forming objects. We therefore classified them using a scheme in- 
tended to differentiate between HERG and LERG radio galaxies. 
Our preliminary classification was visual, i.e. objects were classed 
as Ae' (corresponding to HERGs) if they showed strong high- 
excitation lines such as [Olll], [Nil], [Mgll], [Cm], [Civ] or Lya, 
and as 'Aa' (corresponding to LERG) otherwise, usin g a similar 
classification scheme to that of Mauch & Sadler l2007h - see their 
Section 2.5 for more discussion of this approach and its reliability. 
However, we then found that the equivalent width of the [Olll] line 
gave a very similar division between objects with the advantage of 
removing the subjective element of the visual classification. In the 
final analysis we classified all galaxies with SNR([OlIl])> 3 and 
EW([Om])> 5A as Ae' (HERG-like) and all objects not otherwise 
classified as Aa' (LERG-like). The choice of 5A as the equivalent- 
width cut gives the best match to our preliminary visual analysis, 
but we verified that small variations in this choice made little or no 
difference to the results presented in the rest of the paper. 

A summary of the classification scheme and the number of 
objects in the sample in each of the emission-line classes is given 
in Table Q] Throughout the rest of the paper, we retain a dis- 
tinction between the observational classifications (SF, Aa, Ae, 
AeB) and the physical distinction between star-forming non-AGN 
sources, LERGs and HERGs; we discuss how well the observa- 



3 http://www.mpa-garching.mpg.de/SDSS/DR7/ 



tional emission-line classifications map on to the physical distinc- 
tions in the course of the paper, with a summary in Section|4] 

2.3 Herschel flux-density measurements 

The classification over the GAMA fields and the removal of stars 
gives us 2066 objects, all of which have positions, SDSS identifi- 
cations, FIRST flux densities, spectroscopic redshifts from SDSS, 
GAMA proper or the spare-fibre programme, and spectroscopic 
classifications. Our next step was to extract flux densities for these 
objects from the H-ATLAS 'Phase 1' images. 'Phase 1' of H- 
ATLAS consists of observations of 161 square degrees of the sky 
coincident with the GAMA fields, including the much smaller SDP 
field discussed by H10; further information on the Phase 1 dataset 
will be provided by Hoyos et al. and Valiante et al. (in prep). We 
discarded all GAMA objects which were outside the area covered 
by H-ATLAS (i.e. where flux densities were not available): this re- 
duced the sample to 1836 objects, and it is this 'H-ATLAS subsam- 
ple' that we discuss from now on. 

H-ATLAS maps the FIR sky with HerscheVs Spectra l 
and Photometric Imaging Receiver (SPIRE; iGriffin et all |2010|) 
and the Photodetecto r Array Camera and Spectrometer (PACS; 
IPoglitsch et all l20ld) . The p rocess of deriving th e im ages used 
in this paper is described by IPascale et alj J201 lh and llbar et all 
J2010h for SPIRE and PACS respectively. For each of the objects in 
our H-ATLAS subsample we derived the maximum-likelihood es- 
timate of the flux density at the object position in the three SPIRE 
bands (250 /tm, 350 /im and 500 /im) by measuring the flux den- 
sity from the PSF-convolved H-ATLAS images as in H10, together 
with the error on the fluxes. We also extracted PACS flux densities 
and corresponding errors from the images at 100 and 160 flm us- 
ing circular apertures appropriate for the PACS beam (respectively 
15.0 and 22.5 arcsec) and using the appropriate aperture correc- 
tions, which take account of whether any pixels have been masked. 
We add an estimated absolute flux calibration uncertainty of 10 per 
cent (PACS) and 7 per cent (SPIRE) in quadrature to the errors 
measured from the maps for the purposes of fitting and stacking, 
as recommended in H-ATLAS documentation, but this uncertainty 
is not included when considering whether individual sources are 
detected. 

Only 368 of the H-ATLAS subsample (20 per cent) are de- 
tected in the co nservative '5 a' H-A TLAS source catalogue (created 
as described bv lRigbv et all [2oTlh . This is a similar 5cx detection 
fraction to that obtained by H10. We can relax this criterion for de- 
tection slightly, as we know that there are objects (the host galaxies 
of the radio sources) at the positions of interest. A detection cri- 
terion of 2t7 implies that 2.3 per cent of 'detected' sources will 
be spurious, which is acceptable for our purposes. However, care 
needs to be taken when applying such a criterion to the H-ATLAS 
data. The images at 250, 350 and 500 are badly affected by 
source confusion, and this means that the statistics of the 'noise' 
- including confusing sources - are not Gaussian. We have there- 
fore conservatively determined our 2(7 cutoff by sampling a large 
number of random background-subtracted flux densities from the 
PSF-convolved maps, and determining the flux level below which 
97.7 per cent of the random fluxes lie, to get a flux density limit 
which takes account of confusion. This process returns twice the 
local r.m.s. noise if the noise is Gaussian, which turns out to be the 
case for the PACS data, but gives substantially higher flux density 
limits of 24.6, 26.5 and 25.6 mjy for the 250, 350 and 500-Jlm 
SPIRE maps respectively, corresponding to around 3.8 times the 
local noise estimates for 250 ^m. These limits are essentially inde- 
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Table 1. The classification scheme used in this paper and the number of objects in each class in the H-ATLAS subsample. Also shown are the numbers of 
objects after the application of the 'SF cut' based on the radio/FIR relation, as described in the text. 



Name 


Characteristics 


RL AGN class 


Number in sample 
(Total) ('SF cut' applied) 


Aa 


AGN spectra with EW([OlII])< 5A 


LERG 


1247 1186 


Ae 


AGN spectra with EW([OlII])> 5A 


HERG/NLRG 


199 156 


AeB 


AGN spectra with strong broad high-excitation lines 


HERG/BLRG/QSO 


187 194 


SF 


Star forming galaxy based on BPT or //a-radio correlation 




191 8 



pendent of the local noise estimates (from the noise maps), which 
is as expected since the upper tail of the flux density distribution 
in the maps is dominated by the effects of confusing sources. In 
what follows, we say that a source is 'detected' in a given band if it 
lies above these confusion limits (for the SPIRE data) or above the 
standard 2(7 value (for PACS). By these criteria, 486 sources (26 
per cent) are detected at 250 /im, the most sensitive SPIRE band; 
the number falls to 244 (13 per cent) at 500 flm and 328 (18 per 
cent) at 100 flm. 

We compared this radio-galaxy sample to the sample of VI 2, 
which uses the method described in H10 to select candidate ra- 
dio galaxies, requi ring a cross-match b etween the NVSS and 
the UKIDSS-LAS jLawrence etafl, |2007|) . over the original 135- 
square-degree Phase 1 field. 786 of the current sample match ob- 
jects in the sample of VI 2, and for those objects we find good 
agreement between the NVSS and FIRST flux densities, suggesting 
that there is little missing flux. The objects that are in the H-ATLAS 
subsample but are not identified as radio galaxies in the sample of 
V12 are either not LAS sources or are faint radio sources that fall 
below the NVSS flux density limit but are detectable with FIRST, 
and so would not be expected to be in our NVSS catalogue. We 
conclude that there is good consistency between the method used 
here and the method of H10, in the set of objects where they over- 
lap, and that there is no reason to suppose that the results are less 
robust for the population of faint radio sources that we study for the 
first time in this paper. 

As noted above, our spectroscopically identified sample is not 
complete, in the sense that not all objects that would meet the se- 
lection criteria for spectroscopy have high-quality spectra, and this 
should be borne in mind in what follows. No selection bias has 
been consciously imposed by our choice of objects for spectro- 
scopic analysis. 

2.4 Luminosity and dust mass calculations 

The rest-frame 1 ,4-GHz radio luminosity of the sample sources is 
calculated from the FIRST 1.4-GHz flux density and the spectro- 
scopic redshift, assuming a = 0.8 as in H10. (We comment on con- 
straints on the spectral index of objects in the sample in the next 
subsection.) 

H10 used integrated FIR luminosities, but these depend very 
strongly on the assumptions made about the underlying spectrum, 
in particular the jS and temperature of the modified blackbody 
model which is assumed to describe the data. In this paper we 
instead use the monochromatic luminosity at rest-frame 250 flm, 
L250 This has the advantage that the assumptions we make about 
the spectrum only affect the A"-correction, and so have negligible 
effect at low redshift. We still have to make a choice of the spectrum 
to use for A"-correction, since we cannot fit models to the vast ma- 
jority of our objects. H10 used a modified blackbody with T = 26 
K, j8 = 1.5, but in this paper we use T = 20 K, j3 = 1.8, for reasons 
that will be justified by temperature fits in Section [3~5l 



The disadvantage of this approach is that we lose the ability to 
estimate the star-formation rate directly from the integrated FIR lu- 
minosity, as we attempted to do in H10: however , the relationships 
commonly used to do this (e.g., those given bv lKennicutlll998l) 
are calibrated using starburst galaxies and are not necessarily ap- 
plicable in the temperature and luminosity range that most radio 
galaxies occupy. Instead, we can consider the 250-^ m luminosity 
as representing a dust mass (as in iDunne et alll201 lh : the 'isother- 
mal' dust mass, i.e. the mass derived on the assumption of a single 
temperature for the dust, is given by 

M iso = ^0 (!) 

47CK250B(V250,T) 

where K"25o 1S the dust mass absorption coefficient, which Dunne et 
al. take to be 0.89 m 2 kg -1 , and B(v, T) is the Planck function. It is 
clear for this mass estimation method, and also turns out to be the 
case for the more complex method discussed by Dunne et al., that 
for a roughly constant T we have a linear relationship between mass 
and luminosity, while we also expect a strong correlation between 
L250 and T for a fixed dust mass. Moreover, we expect high values 
of T to be indicators of strong star formation, independent of M[ so . 
In this paper we will initially use L250 to indicate possible differ- 
ences in star formation, and use comparisons of fitted temperatures 
T to confirm them. Later we will show that L250 can be calibrated to 
give a quantitative measure of star-formation rate, subject to some 
important caveats. 

It is important to note that the Herschel SPIRE PSF has a 
FWHM of 1 8 arcsec at 250 flm, which corresponds to linear sizes 
up to ~ 150 kpc at the redshift of the most distant objects in our 
sample. As we noted in H10, the luminosities we measure, and any 
corresponding dust masses or temperatures, apply not just to the 
host galaxy of the radio source but also to its immediate environ- 
ment. Star formation associated with a given AGN might actually 
be taking place in a merging system or a nearby companion galaxy. 



3 RESULTS 

3.1 Subsample properties 

Table \T\ gives the numbers of objects in the H-ATLAS subsam- 
ple that fall into the various emission-line classes defined above. 
We see that absorption-line only or weak emission-line spectra 
('Aa': unambiguously corresponding to the expected spectra of 
'low-excitation' radio galaxies or LERGs) dominate the popu- 
lation. There are then roughly equal numbers of the 'Ae' ob- 
jects, corresponding to the high-excitation narrow-line radio galax- 
ies (HERGs, or NLRGs), broad-line objects ('AeB') and objects 
classed as star-forming on the basis of their spectra ('SF'). 

The redshift distributions within the emission-line classes are 
somewhat different. The objects in the Aa and Ae classes have 
very similar redshift distributions, with median redshifts around 
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Figure 1. Radio sources with a 500- jim detection plotted on the diagnostic 
plot of [Lopez-Caniego e t alj 120121) . Colours indicate different emission- 
line cla sses. The dashed line indic ates the threshold in FIR/radio radio 
used by Lopez-Caniego et aljj2012l) ; below this line, synchrotron emission 
might be bright enough to affect the SPIRE bands. 



0.4, as we might expect for bright galaxies drawn from the par- 
ent (SDSS) sample, and maximum redshift ~ 1. We cannot distin- 
guish between the redshift distributions of the Aa and Ae classes 
on a Kolmogorov-Smirnov (KS) test at the 3o" confidence level. 
The SF galaxies have a clearly different distribution, with median 
z ~ 0.08 and maximum z ~ 0.3, suggesting that these are mainly lo- 
cal, fainter galaxies (as expected from the known different luminos- 
ity fu nctions of the AGN and SF populations; see lMauch & Sadleil 
|2007|) . We retain the SF objects in the sample so as not to exclude 
the possibility, at this stage, that some are NLRG with strong star 
formation. The broad-line objects have a much wider redshift dis- 
tribution, with median z ~ 1.3 and maximum z ~ 3.7. These objects 
are clearly mostly quasars that are in the sample due to their bright 
AGN emission. Similarly, if we consider the radio flux density dis- 
tributions, we cannot distinguish between the Aa or Ae classes at 
high confidence with a KS test, but the SF objects have a signif- 
icantly different flux distribution from the Aa and Ae, tending to 
have fainter radio flux densities. 

The AeB objects are systematically very much brighter in the 
radio, suggesting that the combination of radio and optical selec- 
tion for these quasars is picking up strongly beamed objects, and 
this is true even if we consider only the z < 1 subsample of the 
AeB objects. Among other things, this means that we need to be 
alert to the possibility of non-thermal contamination in the Her- 
schel bands. To check this, we cross-matched the objects in our 
sample to the GMRT catalogue of Mauch et al. (in prep.), who have 
imaged the majority of the Phase 1 area at 325 MHz, using a simple 
positional matching algorithm with a maximum offset of 5 arcsec. 
A total of 536/1836 objects have counterparts in the GMRT cata- 
logue; the low matching fraction reflects the incomplete sky cov- 
erage and variable sensitivity of the GMRT survey, as described in 
detail by Mauch et al. Nevertheless, we can look for spectral index 
differences in the matching objects. The number of cross-matches, 
together with the mean spectral index and the values at the 10th and 
90th percentile, are tabulated as a function of emission-line class in 
Tabled 

While the large number of non-detections in the GMRT sur- 
vey means that we cannot carry out a detailed analysis, we note 



first of all that the mean spectral index of detected sources is close 
enough to our previously adopted value of 0.8 that our ^-correction 
in the radio will not be badly in error, and secondly that the mean 
spectral index of the AeB objects is very much flatter than any of 
the other emission-line classes, although there is still clearly a pop- 
ulation of steep-spectrum AeBs. Given that the detected objects are 
likely to be biased, if at all, towards the steep-spectrum end of the 
intrinsic distribution, it seems likely that the AeB objects contain a 
significant number of flat-spectrum quasars. 

We investigated this issue further by conside ring the diagnos- 
tic methods used by lLopez-Caniego etal] |2012f) in searching for 
blazars. They relied on detections at 500 /im, and, as noted above, 
only a small fraction of our sources have 2o" detections at that 
band. We plotted the sources that do on the diagnostic ra dio/FIR 
colour-colour diagram used bv lLopez-Caniego et alj d2012l) , which 
is intended to search for non-thermal contamination in the SPIRE 
bands; the result is shown in Fig. [TJ We see that of the 24 AeB 
objects with 500-^im detections, about half lie in the region occu- 
pied by the Lopez-Caniego blazar candidates in which synchrotron 
emission might affect the SPIRE bands, a much higher fraction than 
for any other emission-line class. While red 500/350-^im colours 
may just be an indication of low dust temperatures, and the AeB 
sources have higher redshifts than the comparison objects, this is 
a further sign that the AeBs cannot safely be merged with the Ae 
objects in what follows. 

3.2 Herschel and radio luminosity 

Fig. |2]shows the IR luminosity, L250, against the radio luminosity 
for all the objects in the sample. This plot shows several important 
features of the sample. First, we note that the vast majority of the 
broad-line objects (in red) lie at the very high-luminosity end of the 
plot, presumably due to their high redshifts. As we noted above that 
some of these objects may well have FIR fluxes contaminated by 
non-thermal emission, and as their high redshift makes it difficult 
to compare them with radio galaxies in any case, we exclude them 
from further analysis. 

Fig.|2]also shows the expected linear radio-FIR correlation for 
star-forming objects (magenta lines), together with the dispersion 
seen in that relationship, b ased on the parame ter 5250 > which is de- 
fined as logio(L25o/£-i.4) jjarvis et all l201ol) . We see that objects 
classed as SF on the basis of their emission-line properties or their 
radio/Ha relation (black points) almost all lie in this region of the 
plot and close to the best-fitting line; there is some positive devia- 
tion above the line at low luminosities/redshifts, but this was also 
seen by larvis et al.. However, at higher luminosities, a number of 
objects of other emission-line classes also fall in the star-forming 
region, meaning that their radio emission is not bright enough to 
definitively classify them as radio galaxies. Conservatively, every 
object that lies in the star-forming region of this plot should be 
excluded from a discussion of the FIR properties of radio galax- 
ies; following VI 2, we adopt a cut at 17250 > 1-3. The numbers of 
sources remaining, if these objects are excluded, are given in Table 
[T] The vast majority of the SF objects are removed by the cut (here- 
after the 'SF cut'), and although some of the 8 remaining sources 
may be radio galaxies which have erroneously been classified as SF, 
we conservatively exclude them from subsequent analysis (given 
the small numbers involved, including them as though they were 
Ae objects would have little effect on our results). 

Considering only the remaining objects, which we expect to 
be radio-loud AGN, we see that these span a very large range in ra- 
dio luminosity, from 10 22 to 10 27 W Hz -1 if we ignore the broad- 
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Table 2. Spectral indices between 325 MHz and 1 .4 GHz for sample sources with GMRT survey counterparts 



Source type 


Number of matches 


Mean spectral index 


10th percentile 


90th percentile 


All 


536 


0.70 


0.27 


1.09 


SF 


45 


0.88 


0.58 


1.34 


Aa 


356 


0.77 


0.37 


1.09 


Ae 


55 


0.68 


0.34 


0.99 


AeB 


80 


0.50 


-0.03 


1.12 




Figure 2. 250- Jim luminosity against radio luminosity for all the objects in the sample. Stars indicate Herschel detections at 2o" or better as defined in the text, 
points show 2o~ upper limits in IR luminosity derived from the confusion limit. Colours correspond to emission-line classes as follows: Aa, blue; Ae, green; 
AeB, red; SF, black. The solid magenta li nes indicate the ex pected radio-FIR correlation for star-forming objects, <j250 = 1.78, and the approximate scatter 
about this relation, 1 .4 < g 2 50 < 2.1 (from ljarvis et al.U2010t> . The solid orange line shows our adopted 'SF cut' at ^250 = 1-3, and the dot-dashed vertical line 
shows the nominal FRI/FRII break luminosity. 



line objects. The vast m ajority of these lie below the nominal FRI- 
FRII luminosity divide jFanaroff & Rilevi 1 191 4) of 1.2 x 10 25 W 
Hz -1 (plotted on Fig. [2] for reference) and so would normally be 
classed as low-luminosity radio galaxies, though we emphasise that 
the FRI/FRII division is a morphological one and we have made 
no attempt to classify these objects morphologically. In terms of 
our observational emission-line classifications, we see that Aa ob- 
jects dominate numerically by a large factor, but that there are Ae 
objects at all powers. Assuming that Aas trace LERGs and Aes 
HERGs, this is consistent both with what i s seen in brighter radio- 



selected samples at low redshift (see, e.g^ Hardcastle et~all 120091) 



and with the work of iBest & Heckmar] d2012l) over a comparable 
luminosity range. The Aa and Ae objects left after the SF cut has 
been made have redshift and radio luminosity distributions that are 
indistinguishable on a KS test (Fig. |3), but this is not surprising, 
since the differences in the slope of the luminosity function for the 
two populations, leading to the dominance of HERGs at high lu- 



minosities, start to becom e significant only at L\ 4 > 10 25 W Hz 1 
jBest & HeckmaiJ,l2012h . where we have relatively few sources. 



3.3 LERG/HERG comparisons and stacking 

Some differences between the FIR properties of the Aa and Ae ob- 
jects after the SF cut are immediately obvious on inspection of the 
data. For example, 53/156 (34 per cent) of the Ae objects are de- 
tected at the 2a level or better at 250 [lm, while only 93/1186 (8 
per cent) of the Aa objects are detected at this level. 

Overall, both individual sub-samples still being considered 
(i.e. Aa and Ae after SF cut) are significantly detected with respect 
to the background at all three Herschel-SPIRE bands. We follow 
H10 in testing this with a KS test on the distribution of flux den- 
sities compared to random flux densities from the field; the high- 
est null hypothesis probability is 6.6 x 10~ 8 for the Aa sources at 
500 flm, corresponding to a minimum significance of 5.5o" for both 
classes and all SPIRE bands, and the significance is much higher at 
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250 flm. For the two PACS bands, the Aa sub-sample is detected 
at around 98 per cent confidence (i.e. a marginal detection, 2.3c) 
but the Ae objects are significantly detected with a null hypothesis 
probability around 2 x 1(T 7 (5.2CT). 

We are therefore able to adopt the approach of H10 and divide 
our sources into luminosity and redshift bins for a stacking analysis. 
Since we have relatively few Ae sources, we use only three bins in 
both, ensuring that the highest-luminosity bin includes all sources 
above the nominal FRI/II luminosity boundary at L1.4 = 1.2 x 10 25 
W Hz -1 , which, as noted above, can roughly be taken to sepa- 
rate 'low-power' and 'high-lower' radio galaxies. We then used KS 
tests to see whether these subsamples were detected (distinguished 
from the background flux density distribution) at each of the five 
H-ATLAS wavelengths. The results are shown in Tables[3]and[4] 

As found by H10, the detection of all our subsamples is best 
at 250 flm, although with this larger sample most bins are signif- 
icantly detected at 350 flm as well (see Tables [5] and |4j. 500-/im 
detections are less robust, and only the low-redshift Ae subsamples 
are significantly detected in the PACS bands. We therefore rely on 
the 250 flm flux densities for our first estimate of luminosities. With 
the division of the samples into the two emission-line classes, we 
can see that the mean 250-flm flux density for the Ae objects is 
much higher than for the Aa objects in every bin. 

Our stacking analysis follows the method of H10; we de- 
termine the luminosity for each source from the background- 
subtracted flux density, even if negative, on the grounds that this 
is the maximum-likelihood estimator of the true luminosity, and 
take the weighted mean within a bin to estimate stacked bin lumi- 
nosities. Unlike H10, we determine errors on the bins by bootstrap- 
ping, having verified that this method gives very similar results to 
the much more time-consuming and complex method used in the 
earlier paper. As Fig. [4] shows, we find a clear and significant dif- 
ference between the FIR luminosities of the Ae and Aa objects in 
every bin in either radio luminosity and redshift. 

3.4 Comparison with normal galaxies 

As a result of our procedure for generating our radio galaxy cat- 
alogue in H10, we automatically had a comparison galaxy pop- 
ulation. There is no equivalent in the present work, in the sense 
that there is no galaxy population selected in the same way as the 
radio-loud objects. Spectroscopic redshifts for the GAMA sample 



run out at around z ~ 0.6 be cause of the magnitude limit used by 
GAMA jDriver et all l201lh while our spectroscopic sample ex- 
tends to fainter galaxies and higher redshifts. However, a rough 
comparison with radio-quiet galaxies is useful to put our results 
in context. We therefore constructed a comparison galaxy sample 
as follows: 

(i) We based the sample on the galaxy catalogue over the Phase 
1 fields provided as part of the H-ATLAS data release, constructed 
in the manner described bv lSmith et alj J201 lh . and selected galax- 
ies that had either a spectroscopic redshift (from GAMA or SDSS) 
or a photometric redshift with nominal error < 0.1, had measured 
SDSS r and i magnitudes, and were not point-like in r or i. 

(ii) From this sample we took all objects which lay on the ob- 
served H-ATLAS fields and measured their background- subtracted 
250- flm flux densities as for the radio galaxies (giving 318,244 ob- 
jects in total, the vast majority with only photometric redshifts). 
We excluded at this point all objects that formed part of the radio- 
galaxy sample. 

(iii) We A"-corrected the r-band absolute magnitudes of the radio 
gala xies and the comparison sample to z = using KCORRECT v. 
4.2 l lBlanton & Roweisll2007h . 

(iv) Comparing the range of r-band absolute magnitude in the 
radio-galaxy sample with that in the comparison galaxies (Fig. 
we saw that the radio-selected objects tend to be bright galaxies at 
all redshifts. At lower z there is a tendency for t he Ae galaxies to 
be fainter than the Aa (as s een by, for example, iTasse et aT]|2008l 
and lBest & Heckmanll2012l ; note that our sample is not complete, 
which reduces the extent to which we can draw conclusions from 
this observation), but they occupy similarly bright galaxies at high 
z (Fig. IE 

(v) Clearly for even a rough comparison we should compare the 
radio galaxies with optical objects of comparable magnitudes. In 
each of 14 bins of width Az = 0.05 between z = 0.1 and z = 0.8, 
we selected only the comparison galaxies that lay in the absolute r 
magnitude range spanned by the Aa and Ae objects (post-SF cut) 
in that redshift range. 

We then stacked the Herschel FIR luminosities of the galax- 
ies in those 14 bins, deriving them from the 250-/lm flux densities 
on the assumption T = 20 K, /3 = 1.8 in the same way as for the 
radio galaxies. These stacks are plotted as a function of z on Fig. 
[4] We emphasise that this is intentionally a crude comparison: for 
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Figure 4. Comparisons between the IR luminosities of the Aa and Ae objects, after excluding objects near the star-formation line, as a function of (top) radio 
luminosity and (bottom) redshift. Symbols and colours as for Fig. [2] The large crosses indicate the results of stacking the IR luminosities of all Ae (green) and 
Aa (blue) objects in the corresponding radio luminosity or redshift range. Other types of object are plotted but not stacked. The orange bins indicate stacking 
of comparison galaxies, as described in the text. 
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Figure 5. The M^-redshift distribution (X-corrected as described in the text) of sample objects in the Aa (blue) and Ae (green) classes, after the 'SF cut', 
plotted on top of the same information for the comparison galaxies discussed in the text (greyscale shows a density plot with square-root transfer function). 
Red lines show the range of absolute magnitudes of the radio-selected sample, used to generate the comparison galaxy sample discussed in the text. The peak 
in redshift around z ~ 0.8 is probably an artefact of the photometric redshifting, but does not affect our comparison sample. 
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Figure 6. Histogram of (left) the absolute r magnitudes of Aa (blue) and Ae (green) galaxies in the sample, after application of the SF cut, and (right) 
the absolute magnitudes of the comparison objects after the magnitude range selection. Aa objects are generally more massive galaxies than Aes, but the 
distributions have substantial overlap; the distribution of comparison galaxy absolute magnitudes has an intermediate peak. 
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Table 3. Mean bin flux densities and K-S probabilities that the Herschel fluxes of objects in redshift bins (after the SF cut) are drawn from the background 
distribution, as a function of emission-line class and wavelength. Low probabilities (below 1 per cent) imply significant differences between the bin being 
considered and the distribution of flux densities measured from randomly selected positions in the sky, as described in the text. Note that the bins do not include 
quite all the objects in the sample. 



Class 


z range 


Objects 




Mean bin flux density (mjy) 








K-S probability (%) 








in bin 




SPIRE bands 




PACS bands 






SPIRE bands 




PACS bands 








250 jttm 


350 fim 


500 ,ttm 


100 ,tim 


160 


250 ,u 


ii 


350 ^m 


500 llm 


100 fim 


160 llm 


Aa 


0.00 - 0.30 


399 


5.5±0.3 


1.0 ±0.4 


2.5 ±0.4 


2.7 ±1.5 


8.7± 1.9 


< 10" 


3 


< 10" 3 


0.3 


2.2 


3.6 




0.30-0.50 


475 


4.7±0.3 


0.7±0.3 


2.1 ±0.4 


3.1 ±1.4 


7.0±1.7 


< 10" 


-3 


0.004 


0.09 


46.1 


9.9 




0.50-0.90 


310 


5.3±0.4 


1.4 ±0.4 


2.6±0.5 


0.2±1.7 


5.9±2.1 


< 10" 


-3 


< 10" 3 


0.01 


41.1 


35.4 


Ae 


0.00 - 0.30 


62 


27.9 ±0.9 


10.2±0.9 


6.4±1.1 


68.2±3.8 


50.3 ±4.6 


< 10" 


-3 


< 10" 3 


0.05 


< 10 -3 


0.06 




0.30-0.50 


57 


22.8 ±0.8 


8.1±0.9 


3.8±1.1 


25.6±3.8 


30.7±4.7 


< 10" 


-3 


0.005 


2.1 


0.3 


1.5 




0.50 - 0.90 


37 


20.4 ±1.0 


9.8±1.2 


7.6 ±1.4 


23.5±5.0 


14.2±6.1 


< lO- 


-3 


0.002 


0.10 


11.9 


20.8 



Table 4. Mean bin flux densities and K-S probabilities that the Herschel fluxes of objects in luminosity bins (after the SF cut) are drawn from the background 
distribution, as a function of wavelength. Notes as for Table|3] 



Class 


Rang 


e in 


Objects 




Mean bin flux density (mly) 






K-S probability (%) 






it. 4 




in bin 




SPIRE bands 




PACS bands 




SPIRE bands 




PACS bands 










250 llm 


350 ^m 


500 lira 


100 fim 


160 llm 


250 lira 


350 ^m 


500 llm 


100 fim 


160 lira 


Aa 


22.0 


-24.0 


457 


4.7±0.3 


0.5±0.3 


1.9±0.4 


3.2±1.4 


6.9 ±1.8 


< 10~ 3 


0.002 


0.2 


2.1 


14.5 




24.0 


-25.0 


589 


5.7±0.3 


1.5±0.3 


2.6±0.3 


2.4±1.2 


7.5± 1.5 


< 10~ 3 


< 10~ 3 


< 10~ 3 


49.6 


8.1 




25.0 


-28.0 


140 


3.6±0.5 


0.7±0.6 


2.8 ±0.7 


-1.7±2.5 


7.4 ±3.1 


0.006 


0.4 


1.8 


34.0 


19.8 


Ae 


22.0 


-24.0 


71 


19.0±0.8 


6.3±0.9 


3.7 ±1.0 


27.2 ±3.4 


30.7±4.2 


< 10~ 3 


< 10~ 3 


1.7 


0.006 


2.0 




24.0 


-25.0 


57 


35.3 ±0.8 


14.2±0.9 


8.0±1.1 


69.2 ±4.0 


52.5±4.8 


< 10~ 3 


< 10~ 3 


0.01 


0.002 


0.09 




25.0 


-28.0 


28 


15.2± 1.2 


7.2± 1.3 


6.4 ±1.6 


24.3 ±5.6 


7.8 ±7.0 


< 10~ 3 


0.004 


0.4 


4.0 


94.5 



example, we could also have performed a colour selection on the 
comparison galaxies, but this would have involved an investigation 
of the optical colours of the Aa and Ae objects, accounting for pos- 
sible AGN contamination, which we wish to defer to a later paper 
(Ching et al., in prep.); similarly, we are not attempting to sepa- 
rate spirals and ellipticals in the comparison sample, and we have 
made no attempt to match the actual distributions of magnitudes 
(and thus stellar masses) of the comparison galaxies within the ab- 
solute magnitude ranges used (Fig.[6jl. Our optical selection, which 
is required to allow matching to the radio galaxies, also potentially 
biases us against the most strongly star-forming radio-quiet galax- 
ies, which will tend to be more dust-obscured. However, the re- 
sult of the comparison is clear. The luminosities for the comparison 
galaxies tend to lie in between those for the Aa and the Ae radio- 
loud objects; thus, in the redshift range where we have data, Ae 
galaxies are on average more luminous in the FIR than the average 
galaxy of comparable optical magnitude at a given redshift, and Aa 
galaxies are on average less luminous, at least in the low-redshift 
bins where a comparison is possible. We note that the comparison 
sample, though differently selected, is behaving in a very similar 
manner to that of H 10 in terms of its FIR luminosity as a function 
of redshift; but now that we have a larger sample and can separate 
radio galaxies by emission-line class, we are able to see differences 
between the radio-loud and radio-quiet populations. 



3.5 Individual dust temperatures and masses 

Our large sample and the availability of the PACS data allow us to 
investigate the temperatures and isothermal dust masses (Section 
|2.4t of radio-selected objects for the first time. 

We began by fitting single-temperature modified black-body 
models to all the sources for which this was possible. We se- 
lected all objects which had a 2a detection, as defined above, in 
at least two Herschel bands, in order to ensure that there was at 
least in principle sufficient signal-to-noise to constrain parameters, 



and then used standard Levenberg-Marquardt ■/} minimization to 
find the best-fitting values of temperature and normalization for the 
modified blackbody model to the flux densities measured at all five 
H-ATLAS bands. 

One decision that has to be taken here co ncerns the emissivit y 
parameter /3. Earlier work, including H10 and lDunne et alj 1 201 1 ). 
takes this to be 1.5, but work on local galaxies (e.g. lDavies et all 
120121) has obtained good fits (using much better data than available 
to us) with P = 2. Fitting for /3 across the whole sample (marginal- 
izing over temperature and normalization for each object) we find 
that the best fits are found with j8 = 1.8 (the validity of this ap- 
proach will be discussed by Smith et al., in prep). We fixed /3 to this 
value and then fitted for temperature and normalization, determin- 
ing errors by mapping the A% 2 = 2.3 error ellipse (corresponding to 
lo" for two parameters of interest). For these final fits, only individ- 
ual fits with an acceptable x 2 value (defined as a reduced %~ < 2) 
and a well-constrained temperature (AT/T < 2) are considered in 
what follows. This process gives us 385 measured temperatures and 
normalizations, including 128 Aa objects (10 per cent of the total), 
70 Aes (35 per cent) and 170 SF objects (88 per cent); the 'SF 
cut' was not applied to the parent sample. Integrated FIR luminos- 
ity (Lir, integrated between 8 and 1000 microns) and isothermal 
dust masses were then calculated from the fitted temperature and 
normalization. 

Clearly the quantities we measure in this way are expected to 
be biased towards the brightest and hottest objects, but it is still 
instructive to see how they relate to our emission-line classifica- 
tions. Fig. |7] shows the temperature-luminosity and temperature- 
dust mass plots for these objects broken down by emission-line 
class. We exclude the broad-line objects, because of concerns noted 
above about contamination of the FIR bands by non-thermal emis- 
sion. We see what appears to be a bimodal distribution of temper- 
atures, with one set of objects, here seen to be mostly Aas, having 
temperatures in the range 10 < T < 20 K, while the other, com- 
prising most of the SF and Ae objects together with a significant 
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Figure 7. Left: a temperature-luminosity plot for the SF, Aa and Ae sources with individually determined temperatures. Right: stacked, normalized posterior 
distribution functions (PDFs) for temperature fits to the SF, Ae and Aa emission-line classes. 



minority of Aas, has 25 < T < 40 K. The typical error bar on fitted 
temperature (not plotted for clarity) is of order 10 per cent. These 
temperatures gener ally seem real i stic: t he isothermal dust temper- 
atures measured by iDunne et alj d201ll) span the range 10-50 K, 
and, unsurprisingly, our temperatures also lie in this range, while 
temperatures measured for early-ty pe galaxies in the Herschel Ref- 
erence Survey JSmith et alll2012bl) are ~ 24 K, which is compara- 
ble to what we see for the Ae sample. What is more surprising is 
the cold temperatures found for a number of the Aa objects. Inte- 
grated temperatures of the order of 10 K are perhaps just plausible 
for dust in thermal equilibrium with the old stellar population of 
elliptical galaxies, if the dis tribution of dust is chosen carefully; 
iGoudfrooii & de Jond ( fl995h show that temperatures of tens of K 
are expected for dust in the centres of elliptical galaxies, but the ef- 
fective temperature should be an emission-weighted average over 
the dust distribution throughout the galaxy, and the energy density 
in photons falls off very rapidly with distance from the centre of the 
galaxy, so the integrated temperature depends on dust distribution, 
but could be substantially lower than the peak value. However, the 
lowest temperatures found by the fits would imply dust masses of 
up to 10 I0 M ( t,, which is probably not realistic. Inspection of the im- 
ages for the sources with very low temperatures suggests that sev- 
eral of them are the result of flux confusion: given our 2a flux cuts, 
up to 30 of the Aa sources that we have fitted could be spurious 
detections, and while we do not expect the numbers to be this high 
in practice, confusion seems likely to account for a number of the 
sources with the lowest fitted temperatures and highest dust masses. 
In addition, synchrotron contamination of the SPIRE ban ds, which 
is observed in nearby LERGs like M87 l lBaes et allboid) could be 
affecting a few Aas - several fall below the dividing line in the plot 
of Fig.[T] although not all of these will have the flat integrated radio 
spectrum required for non-thermal emission to appear at 500 jlm. 

Another way of investigating the temperature differences be- 
tween the samples, which does not put so much weight on individ- 
ual objects, is to consider the stacked posterior probability distribu- 
tion over T, marginalizing over normalization, for the fitted objects, 
and this is shown in the right-hand panel of Fig. [7] Here a prior of 



5 < T < 55 K is used. We see that this plot reproduces the broad 
trends seen in the left-hand panel: the SF objects have a fairly well- 
defined peak in temperature at around 26 K, the Aes have a broader 
peak at around 30 K, and the Aas span a range between around 10 
and 30 K. Even taking into account possible contamination by con- 
fused sources and/or synchrotron emission in a few cases, it does 
not seem likely that the Aa and Ae sources have the same intrinsic 
temperature distribution. In addition, it is hard to see how this dif- 
ference in the PDFs could be explained by, for example, different 
j6 values for different populations. 

It is natural to interpret the wide range of temperatures in 
the temperature-luminosity plots in terms of two populations of 
dust: (1), a cold dust component which is always present, which 
is essentially in thermal equilibrium with the old stellar population 
(T ~ 15 K might be reasonable for this as an average over the dust 
properties of an elliptical galaxy, as noted above) and whose mass 
scales with the total galaxy mass, perhaps with some redshift de- 
pendence; and (2), a warmer dust component with T ~ 30 K which 
traces current star formation and whose mass and l uminosity are 
primarily an indicator of the star-formation rate (cf . IDunne et all 
l201ll ; ISmithetal.Ll2012ah . In fact, such a two-temperature model 
might help to explain the lack of objects with T ~ 20 K in Fig. [7] 
although there should be objects where the two components con- 
tribute roughly equally to the Herschel SED, these would tend to 
be poorly fitted with a single-temperature model and would be re- 
jected by our fitting procedure. 

The type of bro ad-band fitting to the FIR through ultraviolet 
SEDs carried out by ISmith et all J2012ah is beyond the scope of 
this paper, but we investigated a two-temperature model by fitting 
the same dataset for the normalizations of two modified-blackbody 
models with j8 = 1.8 and fixed temperatures of 15 K and 30 K (nor- 
malizations here are explicitly constrained to be positive to avoid 
trading off negative flux in one component against positive flux in 
another). These models in general fit less well and therefore give 
fewer sources with acceptable fits, presumably because there are 
objects with well-determined temperatures for the warm compo- 
nent that are significantly different from 30 K, but we do note that 
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Figure 8. The relationships between (left) the integrated IR luminosity and (right) the dust mass for objects with fits of both a single-temperature and a 
two-temperature model. 



they provide good fits to a population of objects that are rejected by 
the % 2 criterion for the single-temperature fits and that generally 
have non-zero contributions to the dust luminosity from both hot 
and cold components. Moreover, and importantly for what follows, 
we find a good correlation between the total luminosities for objects 
where these can be obtained using both methods (Fig. [8}, while 
the estimated dust masses show some systematic differences. This 
suggests that, at least for this sample, the total IR luminosity can 
be used without worrying too much about more complex models, 
while the dust mass must be interpreted with a little more care. If 
we interpret the mass of warm dust (or, equivalently for this model, 
its luminosity) in these fits as tracing star-formation rate, then, us- 
ing the results of JSmith et"alll2012a) . the typical SF object in our 
sample has a star-formation rate of ~ 10M P yr _1 , while the most 
luminous Ae objects might have star-formation rates more than ten 
times higher. More detailed methods for estimating star-formation 
rates are discussed in the following subection. 

Finally, we note that the weighted mean of the best-fitting sin- 
gle temperatures of the Aas and Aes, for /3 = 1.8, is 20.3 K. This 
justifies the assumptions we used for A"-corrections in the luminos- 
ity stacking of Section [3~3l 

3.6 L250 as a star-formation rate indicator; comparing 
emission-line classes 

As we noted in the previous subsection, contributions to L250 are 
made both by cold dust (driven by the old stellar population) and 
warm dust (driven by star formation). It follows that neither L250 
nor the integrated Ljr are reliable indicators of star formation rates 
(SFR) in general. However, they should both be usable to estimate 
SFR for an object whose FIR emission is dominated by emission 
from warm dust: these will be the objects whose best-fitting tem- 
peratures are ~ 25 K or more. For objects with a contribution from 
cold dust, the SFR estimated from L250 or Lir will be an upper 
limit. 

To use the quantity that we have used for stacking, Z.250> 
in this way we need to calibrate the relationship between it and 



SFR. We choose to do this by considering the objects classed op- 
tically as 'SF', as, where temperature information is available, all 
of these have FIR temperatures consistent with being dominated 
by star formation (Section |3. 5b . For SF objects with SDSS spec- 
tra estimated star forma tion rates, derived using the methods of 
iBrinchmann et al.l d2004l) (i.e. by model-fitting to the optical emis- 
sion lines and stellar continuum), are available in the MPA-JHU 
database. Cross-matching our SF objects against this database gives 
a sample of 158 objects with both SFR and L250 estimates; we use 
the median likelihood estimates given in the MPA-JHU database, 
as being the most robust, and take half the difference between the 
16th and 84th percentiles as an estimate of the error on the SFR. As 
these objects are all at low redshifts (z < 0.24, and median z = 0.08) 
we need not be concerned about the X'-correction used to derive 
Z,25o- When we plot L250 against SFR derived in this way, we see a 
good correlation (Fig. [9} with a slope that is, by eye, close to unity. 
A Markov-Chain Monte Carlo regression, taking the errors on both 
SFR and L250 into account and incorporating an int rinsic dispersion 
in the manner described bv lHardcastle et alJ j2009T) . gives Bayesian 
estimates of the slope and intercept of the correlation: 

Iogio(i25o/W Hz" 1 ) = 23.64 + 0.961og 10 (SFR/M@ yr" 1 ) 

Although a slope of unity is not ruled out, we will use this slightly 
non-linear relationship in what follows. We emphasise again that it 
is only valid for objects whose FIR emission is known to be domi- 
nated by warm dust heated by star formation. 

As a sanity check on this approach, we can also estimate the 
relationship between SFR and integrated Ltr by using our temper- 
ature fits from Section [TBI The vast majority (143) of SF objects 
with SFR estimates also have estimates of 7ir and thus LrR, and this 
quantity also correlates well with SFR (Fig. [9). Regression gives a 
linear relation 

log 10 (L IR /L Q ) = 9.90+ 1.00iog 10 (SFR/jtf Q yr" 1 ) 

whose normaliz ation is only a fac tor 1.4 away from the standard 
relation given bv lKennicuttl dl998f) , derived for starbursts. Thus we 
can use our IR-derived SFR with reasonable confidence. 
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Applying the L2S0/SFR relation to Fig. [4] we can see that 
the most luminous detected radio galaxies, at around 9 x 10 25 W 
Hz -1 , should correspond to star-formation rates around 250 Mq 
yr , which does not seem unreasonable - these would be radio 
galaxies associated with starbursts - although it should be noted 
that there is a non-negligible uncertainty associated with the L250 
values of these luminous, high-z sources because of the poorly 
known A*-correction. Individual powerful, high-z radio galaxies 
have been associated with star formation at lev els even higher than 
this tearthel et all 1 20121 : 1 Seymour et alll2012h . The mean L250 of 
the most radio-luminous Aes, with Z4.4 > 10 25 W Hz -1 , corre- 
sponds to 15Mq yr , which is well above the SFRs expected for 
normal ellipticals in the local universe. The factor ~ 4 between the 
stacked L250 values for Aes and Aas means that the mean star for- 
mation rate in the latter is at least a factor 4 below that in the Aes; 
'at least' because the temperature measurements suggest that the 
emission from some, and perhaps most, Aas is dominated by cold 
dust. 

3.7 Stacked dust temperatures, masses and SFR 

As noted above, direct estimation of dust temperatures can only 
be carried out for the brightest (and possibly hottest) objects, and 
so might give misleading results if used in the interpretation of our 
stacking analysis of the whole sample. As an alternative, we can es- 
timate mean temperatures for objects in the sample as follows. We 
bin our objects in redshift or radio luminosity as in the previous two 
sections. For each redshift/luminosity bin, we determine the single 
dust temperature that gives the best x 2 fit to the observed flux den- 
sities of every galaxy in the bin, allowing each galaxy to have a free 
normalization (which may be negative) and taking a fixed j8 = 1.8. 
Errors in this fitted temperature are estimated by finding the range 
that gives A% 2 = 1. We can then use the best-fitting temperature 
and normalizations for all the sources to estimate the 250-/xm lu- 
minosity of the bins, determining error bars by bootstrap as before. 
The results of this process are tabulated in Tables [5] and [6] The %- 
fitting gives acceptable, though not particularly good results, as we 
would expect since, from the analysis of Section [3~5l we know that 
there is a wide range of temperatures in each bin. Nevertheless, we 
can attempt to interpret the results. 

Three points are of interest. Firstly, we note that the luminosi- 
ties we estimate are broadly consistent, within the errors, with the 
luminosities estimated from the stacking analysis of Section [33l 
this gives us confidence that the luminosities from the earlier anal- 
ysis are reasonable and that the assumption of a single tempera- 
ture for the A"-correction does not have a big effect on the inferred 
monochromatic luminosities. The luminosity difference between 
the Aa and Ae spectral classes clearly persists in this analysis. 
Secondly, we see that the temperatures are systematically different 
for the two emission-line classes: Aes have systematically higher 
dust temperatures. Thirdly, we can compute isothermal dust masses 
from eq.[T]using the best-fitting temperature and mean luminosity - 
these are of course a complicated weighted mean of the dust masses 
of all the objects in the bin, but still gives us some information on 
the properties of the galaxies. These mean isothermal dust masses 
are tabulated for each bin in Tables[5]and[6] No very strong differ- 
ence between the dust masses for the emission-line classes is seen 
in these mean masses. It therefore seems plausible that the clear 
observed difference in monochromatic FIR luminosity at 250 flm 
between the populations is driven by a difference in dust tempera- 
ture rather than by dust mass. 

Finally, we can attempt to convert the L250 values from this 



fitting into star formation rates, using the results of Section [3~6l As 
already noted, this gives us upper limits if we have reason to sup- 
pose that some of the FIR emission comes from cold dust unrelated 
to star formation. The results of this conversion, given in the final 
column of Tables [5] and [6] must be treated with caution, therefore. 
Since the mean fitted temperatures of the Aes are > 30 K, their 
SFR estimates may be a reasonable estimate of the true mean SFR 
in these systems; the same is not true of the Aas, and so, again, the 
safest interpretation is to say that the mean SFR in the Aes is of 
order a few tens of solar masses per year, and is at least ~ 0.5 dex 
higher than that in the Aas. 

3.8 Radio source sizes 

V12 noted a strong relationship between the FIR luminosity or tem- 
perature of the objects they studied and the radio source size, in the 
sense that larger objects had systematically lower L250 or T. They 
also showed that this discrepancy was largely driven by the most 
massive objects in their sample. 

It is clearly interesting to ask how this result relates to the ob- 
served differences between emission-line classes. One of us (JSV) 
therefore determined the largest angular size of every object in the 
present sample, taking the sizes used by V12 for objects in common 
between the two samples and otherwise making measurements di- 
rectly from the FIRST images. Where objects were unresolved in 
FIRST, an upper limit was assigned, as described by V12. Scaling 
by the angular size distance, this gives the distribution of source 
physical sizes for the current sample. 

An important caveat in this analysis is that the current sample 
is not complete. This is illustrated by the left panel of Fig.[T0] which 
shows the power/linear-size plot (the 'PD diagram') for the current 
sample together with the equivalent plot for the complete and well- 
studied 3CRR sample. For clarity, upper limits are not marked on 
this plot, but it should be noted that all sources with a physical size 
> 40 kpc are actually resolved. We see that all emission-line classes 
are heavily biased towards smaller physical sizes with respect to 
3CRR: of course, this is not a completely fair comparison, since 
the 3CRR objects are the most luminous objects in the radio sky at 
any given redshift, and we might expect lower-luminosity sources 
to be systematically smaller than the most luminous ones. How- 
ever, there are more subtle signs of bias, such as the fact that there 
are more large AeB sources than there are Aes: this arises princi- 
pally because a broad-line object is more likely to have a bright 
radio core and so to be identified with a galaxy or quasar in our 
original selection (though there is an additional effect due to the 
different redshift distribution of AeBs and Aes). The Aa sources, 
which should have the full range of angles to the line of sight, have 
a length distribution intermediate between the Aes and AeBs, as 
expected. This bias towards compact sources, or sources with com- 
pact cores, is particularly problematic for our sample because of the 
selection from FIRST radio images, which resolve out large-scale 
emission: the sample of V12 will be closer to being complete. 

Having said this, it is still possible to investigate the FIR prop- 
erties as a function of length. To do this we apply the SF cut and 
then stack the FIR luminosities as in Section l3~3l binning by length: 
we use only two length bins and the division is set at 40 kpc to 
ensure that all upper limits are in the correct bin. The results are 
shown in Fig.[H)](right panel). We see first of all that the Ae/Aa dif- 
ference persists in this analysis: Ae sources have higher L250 than 
Aas irrespective of length. Secondly, we see no evidence for any 
length dependence of the L250 of the Ae population, although the 
error bars are large because the sample is small. Thirdly, we note 
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Figu re 9. The relationsh ip between star-formation rate (SFR), determined for the SF class from the MPA-JHU emission-line database using the methods of 
iBrinchmann et al.ll2004l) . and the two types of FIR luminosity discussed in this paper, L250 (left) and Lir (right). The green lines show the results of a MCMC 
regression, as described in the text. Error bars are not plotted for clarity - the uncertainty on SFR can be ±0.5 dex. 



Table S. Results of temperature fitting, dust mass and star formation rate estimation in redshift bins. 



Class Range in Objects Best-fit T Reduced log 10 (L 25 o/W Hz~') log lo (M iso /M ) log lo (SFR/M yr _1 ) 







in bin 


(K) 


x 2 








Aa 


0.00-0.30 


398 


17 +IW 
1 '- u -0.8 


1.65 


23.9 ±0.1 


7.6 ±0.1 


0.2±0.1 




0.30-0.50 


472 


17 4+0.6 
17 -^-0.6 


1.52 


24.3 ±0.1 


8.0 ±0.1 


0.7±0.1 




0.50-0.90 


309 


27.8;" 


1.44 


24.6 ±0.1 


7.7±0.1 


1.0 ±0.1 


Ae 


0.00 - 0.30 


62 


32.8ig 


2.11 


24.3 ±0.1 


7.3 ±0.1 


0.7 ±0.1 




0.30-0.50 


55 


30.9j{2 


1.45 


24.8±0.1 


7.8±().l 


1.2±0.1 




0.50 - 0.90 


36 


34.2l|j 


1.99 


25.1 ±0.1 


8.0±0.1 


1.5 ±0.1 



a marginally significant difference between the L250 values for the 
small and large Aa sources, in the sense noted by V12: the null hy- 
pothesis that these two are equal can be rejected at the 95 per cent 
confidence level. If, instead of stacking Z/550, we fit temperatures 
to all the sources in each bin as described in Section [3771 we find 
no significant difference in luminosities as a function of length for 
either emission-line class, but there are significant differences in 
best-fitting temperature (29.4 ± 0.5 vs 22 ± 2 K for Ae; 21.8 ± 0.5 
vs 12.2 ±0.4 K for Aa), in both cases in the sense that the larger 
objects have lower best-fitting temperatures. This is again broadly 
consistent with the results of VI 2. 



4 DISCUSSION 

The results of the previous section show very clearly that there is 
a difference between the average far-infrared properties of radio 
galaxies whose spectra show strong emission lines and those of 
radio galaxies that do not. How can we interpret this? 

It is first of all important to consider the issue of possible AGN 
contamination, discussed in Section[T] in more detail. There are two 
possible sources for this: (1) emission from the warm dusty torus, 
which is expected to be seen predominantly in HERGs, and (2) syn- 



chrotron emission from the jets and lobes, which may appear in all 
objects. The first of these is particularly important, as torus-related 
emission in the FIR bands might give rise to a HERG/LERG dif- 
ference, but we are confident that it is not a significant effect in our 
sample, for several reasons. Firstly, when the required mid-IR data 
are available, which is not the case for our sample at present, de- 
compositions of the SEDs of radio-loud and radio-quiet AGN tend 
to show that observer-frame Herschel SPIRE bands are dominated 
by cool dust rather than by th e torus component, even for powerful 
AGN with luminous tori (e.g lBarthel et alll2012llDel Moro et all 
120121) . Secondly, we know, as pointed out by H10, that the mid- 
IR torus luminosities even for the most powerful HERGs in our 
sample, if they follow the correl ation between radio po wer and 
mid-IR luminosity established by lHardcastle et alj u2009h . should 
be about two orders of magnitude less than the total FIR luminosi- 
ties estimated e.g. in Fig. [7] implying that even if the torus SEDs 
were strikingly different from those of known radio-loud AGN, it 
would still be energetically impossible for them to affect the ob- 
served FIR emission significantly. The second possible source of 
contamination, synchrotron emission, we believe to affect mostly 
the broad-line objects, as discussed in Section [3~71 together with at 
most a very few of the Aas; there is certainly no reason to expect 
that it would give rise to the observed Aa/Ae difference, since the 
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Table 6. Results of temperature fitting, dust mass and star formation rate estimation in luminosity bins. 



Class 


Range in 
log 10 (Li. 4 ) 


Objects 
in bin 


Best-fit T 
(K) 


Reduced 

X 1 


log IO (L 25 o/W Hz" 1 ) 


log 10 (M iso /M G ) 


log 10 (SFR/M :; yr-') 


Aa 


22.0 - 24.0 


456 


15 \ +0b 

-6.5 


1.50 


24.0 ±0.1 


7.9±0.1 


0.4±0.1 




24.0-25.0 


589 


23.8+2-2 


1.65 


24.5 ±0.1 


7.7 ±0.1 


0.9±0.1 




25.0-28.0 


140 


15-3$ 


1.48 


24.7 ±0.1 


8.6±0.1 


1.1±0.1 


Ae 


22.0 - 24.0 


71 


26.31°;? 


1.70 


24.3 ±0.1 


7.5±0.1 


0.7 ±0.1 




24.0-25.0 


57 


30.9;»J 


2.08 


25.0±0.1 


8.0 ±0.1 


1.5±0.1 




25.0-28.0 


28 


27.0l|] 


1.59 


24.9±0.2 


8.0±0.2 


1.3 ±0.2 
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Figure 10. The physical sizes of radio sources in our sample and their relationship to radio and FIR luminosity. Left: the power/linear-size plot for the sample, 
broken down by emission-line class, and without the SF cut. Colours are as in Fig. [2] for comparison, 3CRR sources are plotted in magenta. Right: the FIR 
luminosity/linear size plot, with stacking in bins at < 40 and > 40 kpc, after the SF cut. Colours and symbols as in Fig. [2] upper limits in length are denoted 
by arrows. 



radio fluxes and spectra of these two classes are very similar. We 
therefore consider it safe to discuss our observations in terms of a 
difference in the properties of cool dust in the two populations. 

Our result cannot be significantly affected by contamination 
by pure star-forming objects whose radio emission is bright enough 
to cause them to be misidentified as radio galaxies. While our spec- 
troscopic classification alone does not identify all objects whose 
radio emission is dominated by star formation (Fig.[2j the combina- 
tion of optical spectroscopy and the 'SF cut' that we impose on the 
radio-FIR luminosity plane, where a clear star-forming sequence is 
visible, should remove all such objects. Moreover, the highest radio 
luminosities in our sample are well above even the~ 10 25 WHz -1 
expected from a starburst of a few thousand Mq yi~ 1 , and we see a 
clear difference between the different emission-line classes in this 
luminosity range. 

Along similar lines, we do not believe that the relationship 
between emission-line class and FIR emission can be a result of 
optical emission-line activity due to the star formation process it- 
self. Our emission-line classification uses [Oin], which, at least at 
high luminosities, is widely used as an AGN indicator, although it 
can be produced by hot, young stars. We do not have a direct esti- 
mate of the [Om] luminosities of our sample objects in the version 
of the GANDALF-derived database we use, but we have derived a 
rough indicator of luminosity from the measured equivalent widths 
and the A"-corrected absolute magnitude in the SDSS g band. Cal- 
ibrating this indicator using the MPA-JHU emission-line measure- 
ments, for which both equivalent width and [Oin] flux are tabu- 



lated, we see that almost all the Ae objects would be expected to 
have L ([Om1) > 10 4Q erg s~*, and are thus in the range classified 
by e.g. iKauffmann et alj ( r2003l) as 'strong AGN'. Further work in 
this area will require measurements of the [Om] fluxes, and ide- 
ally those of other lines, for a large radio-galaxy sample, but we are 
confident that we are assessing genuine AGN activity in the vast 
majority of cases. 

We can therefore move on to interpreting the relationship as 
being one between FIR properties of the host galaxy and the AGN- 
related emission-line properties of radio galaxies that we discussed 
in Section [T] with the Aa population corresponding to LERGs and 
the Ae population to HERGs. It then appears that HERGs, on av- 
erage, have significantly higher L250 than LERGs (Section 13-31 ); 
moreover, HERGs appear to have higher L250 than normal galax- 
ies of comparable absolute magnitude at all redshifts (Section r3.4t . 
In a simple isothermal model, higher luminosity can arise either 
because of higher masses of dust or higher temperatures; what we 
see from the analysis of Sections 13.51 and 13.71 is that it is plausi- 
ble that the dust masses of the different systems are similar, but 
that the mean isothermal temperatures of the HERGs are higher. 
In resolved local galaxies, it has been shown that low-temperature 
dust emission (T ~ 15 K) is driven by the old stellar population, 
while signifi cantly hotter tempe r atures are seen from star-forming 
regions (e.g. iBendo et all |201(]| : iBoquien et all l201lh . By far the 
most obvious interpretation of our result is therefore that the star- 
formation rates are significantly higher in the HERG subsamples 
than in the LERGs, giving rise to a significant component of emis- 
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sion from hot dus t which raises the is othermal temperature as seen 
in the analysis of iDunne et alj d20 1 lh - If so, this is strong confir- 



from the earlier work discussed in Section[T](e.} 
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calibrating L250 as a star formation indicator using SFRs derived 
from local radio-loud star-forming galaxies (Section [3.6l > we have 
been able to quantify this, showing that the mean SFR in the most 
luminous/high-z Aes is probably at the level of around 30Mr., yr _1 , 
and is at least ~ 0.5 dex higher than that in the Aas at all redshifts 
and radio luminosities. 

What does this tell us about the association between star for- 
mation and AGN activity in radio galaxies? The first point to 
note is that the association between a HERG classification and in- 
creased FIR luminosity (and therefore star-formation rate) is statis- 
tical only. It can be seen from Fig.|4]that there are individual LERGs 
with high FIR luminosities, while at the same radio luminosity we 
see HERGs with FIR luminosities 0.5-1 dex lower. Similarly, the 
temperature analysis of Section [331 shows that there are LERGs 
with best-fitting temperatures comparable (within the large errors) 
to those of the warm dust in known star-forming galaxies. Nothing 
appears to require HERGs to be associated with high star-formation 
rates or LERGs to be associated with completely quiescent galax- 
ies, consistent with the conclusions of lTadhunter et alj J201 ll) . This 
suggests that the mechanism of the association is not the simplest 
possible one, in which some single event, such as a merger, always 
triggers both HERG activity and star formation. If this were the 
case, we would not see individual LERGs with high star-formation 
rates (setting aside the possibility, which we regard as remote, that 
these sources are all misidentified HERGs). 

Another piece of evidence supporting this picture comes from 
the lobe length analysis of Section [3~8l If HERGs were associated 
with AGN triggering following a merger, we might expect them 
to show a very strong relationship between FIR properties and lobe 
physical size, since star formation would be expected to peak on av- 
erage at early times in the radio source's lifetime. This was a pos- 
sible interpretation of the results of VI 2, who showed that larger 
sources in general have lower FIR temperatures and luminosities. 
However, our analysis shows that this result is not driven by the 
HERG (Ae) population, and in fact for our sample is more obvious 
for the LERGs (which, however, have considerably better statis- 
tics). 

We would therefore argue that we are not seeing a sim- 
ple triggering relationship, but rather that the difference between 
HERGs and LERGs is that HERGs tend to inhabit environments 
in which star formation is favoured relative to the general galaxy 
population, while by contrast star formation is disfavoured in the 
environments of LERGs. The FIR differences as a function of 
source length would then be explained by some other process, 
such as jet-induced star formation when the bow shock of the 
source is within the host galaxy, which can in principle take place 
in both emission-line classes (although we note there is not yet 
any direct evidence for this process affecting emission seen in the 
FIR band). Such a model is consistent both with all the obser- 



vations to date fe e. iBaldi & Capettil 120081; iHerbert et all l20ld : 
iBest & Heckmanj . |2012| ; Ijanssen et all I20T2I) and with the expla- 



nation of the HERG/LERG dichotomy in terms of accretion mode 
discussed in SectionQ] It will be of great interest to see whether this 
result is confirmed by the larger samples that will be made avail- 
able by the full H-ATLAS dataset, whether it can be extended to 
higher redshifts using deeper spectroscopic or photometric surveys, 
and whether the same results are obtained when the HERG/LERG 



classification is made using data at other wavebands (e.g. X-ray or 
mid-IR). 



5 SUMMARY AND CONCLUSIONS 

The key results from the analysis and discussion above can be sum- 
marized as follows: 

• We have used individual measurements and stacking analy- 
ses to determine the FIR properties (mean luminosities and tem- 
peratures) of a large sample of radio-selected sources with spec- 
troscopic redshifts and HERG/LERG classifications from optical 
spectroscopy. Sources near the known FIR-radio correlation are ex- 
cluded from our analysis; the vast majority of the objects we study 
should be bona fide radio galaxies. 

• We find a clear difference between the FIR properties of the 
two populations in the sense that the rest-frame 250-^im luminosi- 
ties are systematically higher in the HERGs than in the LERGs; 
the host galaxies of LERGs in fact occupy galaxies with lower FIR 
luminosities than normal galaxies matched in absolute magnitude, 
while HERGs tend to have higher FIR luminosities. This difference 
is apparent at all redshifts and all radio luminosities sampled by our 
targets. 

• A comparison of the temperatures and dust masses of HERGs 
and LERGs, stacked in coarse bins, suggests that the dust masses 
are reasonably comparable for the two samples but that the tem- 
peratures in the HERGs are systematically higher. This provides 
strong evidence that the higher FIR luminosities we are seeing im- 
ply, on average, higher star formation rates (which are required to 
raise the mean temperature of the dust) rather than just higher dust 
masses. The low mean temperatures seen for LERGs are consistent 
with what would be expected for quiescent dust which is in thermal 
equilibrium with the photon field of the old stellar population of the 
host galaxy, although the fact that these objects are detected at all 
implies that large masses of dust are present. 

• Quantifying the SFR by calibrating L250 as a star-formation 
indicator in the 'SF' sources known to be dominated by hot dust, 
we find that the mean SFR in the radio-luminous Aes is ~ 30M,;, 
yr , and is at least ~ 0.5 dex higher than that in the Aas at all 
luminosities and redshifts. 

• Consistent with the results of V12, we find that both emission- 
line classes in our sample show some evidence for a dependence of 
FIR properties on radio source size. 

• We argue that there is certainly not a simple triggering re- 
lation, and not even a one-to-one association, between enhanced 
star formation and a particular AGN type (a conclusion con- 
sistent with detailed studies of starburst radio galaxies such as 
that of iTadhunter et aL 1 1201 ]]). However, the statistical trend for 
HERGs to have higher star formation rates is consistent both 
with what is known from o ther wavebands (e.g. IBest & Heckmanl 
120121; Ijanssen et all 120121) and with the general class of models 
1 Hardcastle et all 2007 ) in which HERG activity takes place in 
lower-mass galaxies where the black hole is able to accrete sig- 
nificant quantities of cold gas. 



As more H-ATLAS data and supporting optical imaging and 
spectroscopy become available we expect to extend this work to 
much larger samples, allowing more detailed binning and tem- 
perature analysis, to investigate different methods of carrying out 
the LERG/HERG classification, and to consider radio galaxies at 
higher redshifts in order to search for evidence of cosmological 
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evolution of the star-formation properties of the radio-loud AGN 
population. 
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